Tuning the Ultrafast Spin Dynamics in Carrier-Density-Controlled Ferromagnets 
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Ultrafast strengthening or quenching of the ferromagnetic order of semiconducting Eui_ v Gd, v O was achieved 
by resonant photoexcitation. The modification of the magnetic order is established within 3 ps and detected via 
optical second harmonic generation. A theoretical analysis shows that the response is determined by the inter- 
play of chemically and optically generated carriers in a nonequilibrium scenario beyond the three-temperature 
model. General criteria for the design of spintronics materials with tunable ultrafast spin dynamics are given. 
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The demand for an ever-increasing density and speed of 
manipulation in magnetic information storage has triggered 
an intense search for ways to control the magnetic moment by 
means other than magnetic fields. Following the pioneering 
demonstration of light-induced magnetization change on the 
sub-ps time scale [1], the optical manipulation and control of 
magnetic order by ultrashort laser pulses has developed into 
an exciting research topic in modern magnetism [2-6]. Most 
of the time-resolved investigations are focused on ultrafast de- 
magnetization [1, 7-9]. This involves quenching of magnetic 
order and is often described by a three-temperature model 
(3TM), where the coupled electron, spin and lattice subsys- 
tems are each assumed to be in equilibrium at their respective 
temperatures [1, 7]. However, in special cases, the photoex- 
citation leads to the nonthermal generation and/or enhance- 
ment of magnetic order by photo-controlling the exchange in- 
teraction. For example, the ultrafast generation of ferromag- 
netic order via double-exchange interaction in a strongly cor- 
related manganite [10] and the transient enhancement of fer- 
romagnetism via p-d exchange interaction in a diluted mag- 
netic semiconductor [11] have been observed. These obser- 
vations arouse the question if materials, material parameters, 
and eventually, general criteria can be identified that allow us 
to enhance or attenuate the magnetic order of a system at will. 

In this Letter we show that, depending on the carrier density 
n c , ultrafast strengthening or quenching of the ferromagnetic 
order in Eui_ r Gd T can be achieved via resonant photoexci- 
tation. The change of the magnetic order is established within 
3 ps and detected via magnetization-induced optical second 
harmonic generation (MSHG). The spin dynamics cannot be 
explained by a 3TM, but is a consequence of the photoinduced 
non-equilibrium carrier distribution, as the theoretical analy- 
sis shows. In contrast to transition-metal, rare-earth, or cer- 
tain [11] diluted magnetic semiconductor ferromagnets with a 
high carrier concentration, the sub-ps carrier dynamics in low- 
doped Eui_ v Gd v O is far from equilibrium due to a substan- 
tially longer electronic thermalization time. The resulting spin 



dynamics is explained by the interplay of chemically and op- 
tically generated carriers and the dynamic renormalization of 
the RKKY-like exchange coupling. The change from ultrafast 
enhancement to quenching of ferromagnetism is marked by 
the crossover from the predominantly semiconducting to the 
predominantly metallic doping range of Eui_ x Gd T 0. Aside 
from giving a generalized insight into ultrafast spin as well as 
charge dynamics, a key requirement for spintronics technol- 
ogy, our investigation provides a guide for the selection and 
design of materials for ultrafast optical control of magnetic 
order. 

Gd-doped EuO was chosen as candidate material because 
of the multitude of extreme magnetism-related properties sur- 
passing those of the diluted transition-metal-based magnetic 
semiconductors (e.g., full spin polarization, a simultaneous 
ferromagnetic and insulator-metal transition, colossal magne- 
toresistance, giant magneto-optical effects [12-19]). Undoped 
EuO is a prototype Heisenberg ferromagnet. The magnetic 
moment of 7 jig arises from the strongly localized, half-filled 
4/ shell of the Eu~ + ions on a cubic rock salt structure [inset 
of Fig. 1(a)]. However, because of the strong localization of 
the 4/ orbitals, the direct 4/-4/ exchange coupling with the 
12 nearest Eu 2+ neighbors is too weak to induce ferromag- 
netic order at a Curie temperature as high as Tc = 69 K. In- 
stead, the ferromagnetic order is driven by a virtual exchange 
mechanism [20]: A virtual magnetic exciton is created by a 
localized 4/ electron fluctuating into the empty Eu 5d state so 
that the magnetic exchange coupling J ex can be mediated via 
the spatially more extended 5d orbitals. 

In Eui_ x Gd x O, J ex and Tc are further enhanced by truly 
populating the 5d orbitals with the extra 5d electron intro- 
duced by replacing Eu 2+ by the otherwise magnetically equiv- 
alent Gd 3+ [21, 22], as shown in Fig. 1(a). Note that these 
dopant electrons do not form excitons, but rather occupy the 
Gd 5d impurity orbitals, which in the metallic phase below Tc 
merge with the Eu 5d6s conduction band [23] to produce a 
long-range RKKY contribution to J ex . 
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FIG. 1: (color online) (a) Dependence of Tq (solid triangles) [22] 
and of the square root of the static MSHG intensity V/ MSHG on the 
carrier density n c at T = 4.2 K in Eui_ A Gd x O. Insets show the Gd 
concentration dependence of n c [22] and the crystal structure of EuO. 
(b) Spin-dependent electronic structure of ferromagnetic EuO. The 
Eu 5dbs conduction band is split by the 4f-5d exchange interaction. 
The pump pulse excites the localized 4/ electrons to form 4f-5d 
magnetic exciton states. In the figure the energy splitting of the ex- 
cited 4f b 5dt2 g states and the relevant magnetic exciton states [20] are 
neglected for simplicity, (c) Temperature dependence of V/ MSHG for 
non-doped and 2.65% Gd-doped EuO. 



This mechanism suggests that, alternatively, an ultrafast en- 
hancement of J ex may be driven by resonant optical pump- 
ing of electrons from the Eu 4/ 7 ground state to the 4f 6 5dt2 g 
magnetic exciton state [Fig. 1(b)], since this turns the virtual 
magnetic exciton into a real one, promoting the magnetic or- 
der. Thus, the resonant optical pumping of the magnetic ex- 
citon is an ideal way to control the magnetic coupling in EuO 
on the ultrafast time scale. 

In order to study the ultrafast spin dynamics, epitaxial 
Eu 1 _ v Gd v O(001) films (x = - 19.5%) with a thickness 
of 35 nm were grown adsorption-controlled on two-side- 
polished YAlO3(110) single-crystal substrates by molecular 
beam epitaxy [22]. The Gd concentration x was determined 
by prompt-gamma activation analysis and x-ray absorption 
spectroscopy (XAS) [22]. Uniform growth without secondary 
phases was confirmed by x-ray diffraction. The films were 
protected against air by an amorphous silicon cap layer of 
10 — 20 nm and found to be free of oxygen vacancies within 
the resolution limit of XAS [22]. Due to the strong shape 
anisotropy all Eui_ A Gd x O films possess an in-plane magnetic 
easy axis. Physical properties of our Eui_ A Gd x O films have 
been reported in Ref . [22] . 

The experimental pump-probe setup is shown in the up- 



per inset of Fig. 2. The output of a Ti:sapphire regenerative 
amplifier system operated at 800 nm (E p i, = 1.55 eV) with 
a pulse width of 130 fs and a repetition rate of 1 kHz was 
divided into two beams. One was used for resonant pump- 
ing of the 4/ 7 -> 4f 6 5dt 2g transition at 1.55 eV. The other 
was converted into the frequency-tunable probe beam by an 
optical parametric amplifier. The focus size of the linearly 
polarized pump and probe pulses was about 1 and 0.4 mm, re- 
spectively. MSHG was used for probing the time-evolution of 
the ferromagnetic order. A light wave of frequency co is inci- 
dent onto the sample where it generates a polarization P(2co) 
coupling linearly to the magnetically ordered state. It does 
so with high symmetry selectivity and background-free which 
makes MSHG an ideal probe of magnetic order [18]. The 
MSHG yield / MSHG oc \P(2(o)\ 2 was measured in a normal- 
incidence transmission geometry under an in-plane magnetic 
field of ~50 mT applied in the Voigt configuration for obtain- 
ing a magnetic single-domain state. 

As an initial step we have to clarify the coupling of the 
MSHG signal to the magneti c order. First, as seen in Fig. 1(a), 
the dependence of V/ MSHG on n c follows that on Tq which 
is, in turn, proportional to J ex . (The proportionality between 
Tq and J ex is derived from mean-field theory which has previ- 
ously been shown to describe even the non-Brillouin-like tem- 
perature dependence of the Gd-doped EuO [23].) In the sim- 
plest approximation one may cast all electron dynamics and 
doping dependence into an effective exchange coupling J ex . 
Although this is a non-trivial assumption, it is experimentally 
verified well by Fig. 1(a). We may therefore approximate the 
relation between \// MSHG and J ex (resp. 7c) as linear. Second, 
the temperature dependence of V/ MSHG shown in Fig. 1(c) is 
well reproduced by that of the spontaneous magnetization M 
[18, 22]. Hence, these results lead us to the phenomenological 
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FIG. 2: (color online) Time evolution of the change of the MSHG 
intensity at T = 77 K for EuO doped with 2.65% Gd. The excitation 
density of the pump pulse is ~380 /iJ/cm . The insets show the 
schematic of the time-resolved MSHG spectroscopy and a magnified 
view of the ultrafast increase of the MSHG signal, respectively. Solid 
lines are guides to the eye. 
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Note that in general the effective coupling J ex may be medi- 
ated by the virtual (undoped ground state) or the real (Gd- 
or photo-generated) occupation of the 5d orbitals. Previous 
experiments revealed that at 2h(0 = 2.60 eV the interference 
by temperature-dependent linear absorption effects is avoided 
[18]. We therefore chose it as MSHG probe energy. 

Figure 2 shows a typical time evolution of the photoinduced 
change of the MSHG intensity, normalized by the value be- 
fore optical excitation, A/ MSHG (Af )// MSHG (0). The data were 
taken at 77 K on EuO doped with 2.65% Gd. Following the 
photoexcitation at At = we observe a continuous increase 
of the MSHG intensity up to +30% at At = 3 ps. This is fol- 
lowed by a continuous decrease passing zero at At = 40 ps 
and —25% at At = 3 ns. We conclude that the photoexcita- 
tion populates the 5d orbital, resulting in an enhancement of 
the effective Af-5d exchange and, by aligning the thermally 
fluctuating 4/ spins, of the magnetization M. 

The intra-atomic 4f-5d exchange energy of ~ 0.1 eV [20] 
corresponds to a time of ~ 40 fs so that an instantaneous 
MSHG increase after photoexcitation might be expected in 
contrast to the strikingly different time of 3 ps. As the re- 
cent theoretical analysis shows [24], however, the magnetic 
exchange coupling is mediated via the RKKY interaction. The 
time To required to coherently establish its enhancement in the 
photoexcited region can be estimated from the time it takes the 
electronic correlation to spread through the crystal. The rel- 
evant distance is the RKKY wavelength, and the propagation 
velocity is the group velocity of the magnetic exciton. The cal- 
culated buildup time To of ^10 ps [24] is in surprisingly good 
agreement with the measured value of 3 ps in Fig. 2. The 
ensuing decrease of MSHG intensity on the ns time scale is 
associated with the demagnetization by the heating of the spin 
system through the transfer of the optical excitation energy 
from the electron system to the lattice and subsequent spin- 
lattice relaxation. The observed slow relaxation is consistent 
with the empirical thermal demagnetization time [9] inferred 
from the magnetocrystalline anisotropy constant in EuO [25], 
and the decrease of the MSHG intensity is consistent with the 
laser-pulse heating calculated from the absorption coefficient 
and the heat capacity of EuO [26]. 

As the central part of the present work, we investigated the 
relation between the magnetic coupling dynamics and the car- 
rier density n c in Eui_ v Gd x O In order to simplify this inves- 
tigation, the corresponding dynamics were measured at 10 K, 
where the magnetic moment is saturated. Hence, according 
to Eq. (1), all the non-thermal changes of the MSHG intensity 
can be associated with changes of Jex. Figure 3 shows the tem- 
poral evolution of the normalized change of MSHG intensity 
induced by optical pumping for Gd concentrations between 
and 19.5%. A common behavior of all samples is the decreas- 
ing MSHG yield on a ns time scale. As mentioned above, this 
is caused by the thermal destabilization of the magnetic order 
mediated by the spin-lattice relaxation. In contrast, the MSHG 
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FIG. 3: (color online) Time evolution of the change of the MSHG 
intensity at T = 10 K for Eui_ v Gd r O as a function of time delay 
and Gd concentration. The excitation density of the pump pulse is 
~130 jllj/cm 2 for all traces. The inset shows the magnified view of 
the ultrafast increase and decrease of the MSHG intensity for EuO 
doped with 0.64% and 9.63% Gd. 



response on the ps time scale depends strongly on the Gd 
concentration. At the lowest doping an ultrafast increase of 
the MSHG signal is observed. It strikingly confirms that our 
magneto-optical probe process reflects a variation of J ex rather 
than of the already saturated magnetization M, in contrast to 
the results of earlier time-resolved magneto-optical rotation 
measurements on undoped EuO [27]. By increasing the Gd 
concentration, the magnitude of this initial increase becomes 
progressively more pronounced and reaches a maximum near 
0.25% Gd concentration. However, further Gd doping reduces 
its magnitude, and at high Gd concentration (> 5%) an ultra- 
fast decrease of the MSHG signal is observed. In Fig. 4(a) 
we summarize the pump-induced change of the MSHG yield 
as a function of n c at the fixed delay At = 3 ps. Here, n c has 
been determined for each Gd concentration x from the relation 
displayed in the inset of Fig. 1(a). 

The growth of the MSHG response A/ MSHG // MSHG with 
increasing Gd concentration on the ps time scale is striking 
and surprising at first: According to the widely used 3TM 
[1, 7], the optical pump pulse would, in addition to exciting 
the localized 4/ electrons and creating the 4/ 6 5<if2? magnetic 
excitons, lead to the excitation of the chemically doped car- 
riers in the 5d6s conduction band and to their thermalization 
at the elevated temperature corresponding to the deposited en- 
ergy. This heating of the electron system would invariably 
reduce J ex since the magnetic coupling with electrons in the 
energetically higher band states becomes increasingly antifer- 
romagnetic because of their increased wave number. (This 
was checked by explicitly calculating the RKKY interaction 
for the corresponding electron distribution in the conduction 
band [24].) Thermalization effects would further contribute to 
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FIG. 4: (color online) (a) n c dependence of the photoinduced change 
of the MSHG intensity at At = 3 ps and at T = 10 K extracted from 
Fig. 3. The line is a guide for the eye. (b) Theoretical modelling 
of the RKKY-related photoinduced change of J ex , see text, (c) Sim- 
plified sketches comparing non-thermal and thermal photoexcitation 
dynamics. Left: Photoexcitation of the localized 4/ electrons in low- 
doped Eu A -GcUO. Two-body photon absorption within the 5dbs con- 
duction band is kinematically forbidden. Right: Photoexcitation and 
ultrafast fhermalization of conduction electrons in a metallic ferro- 
magnet, like Gd, due to three-body collisions. Because of the large 
band gap transitions from the 4/ state are energetically inaccessible. 
Blue: Intrinsic or chemically doped, red: optically generated charge 
carriers. For clarity, only the majority (spin-up) earners are shown. 
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with n c [Fig. 4(a)] is that in the low- 
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doped Eui_ x Gd x O the carrier thermalization time x c is 
enhanced over the one in dense itinerant ferromagnets, like 
Ni or Gd, by a factor n^ d /n c sa 100. This is because charge 
thermalization occurs by electron-electron scattering whose 
rate is proportional to the carrier concentration /? Gd in Gd or 
n c in Eui_ x Gd x O. Taking the carrier thermalization time in 
pure Gd as 50—100 fs [7], we thus estimate T c in low-doped 
Eui_ T Gd x O as % c rj 5 — 10 ps. This means that the dynamics 
for the first few ps is dominated by a non-equilibrium charge 
carrier distribution — the 3TM is not applicable. More 
importantly, since electron-electron interactions are scarce 
in this time range, photoexcitation of the chemically doped 
carriers within the 5d6s band into higher band states is 
kinematically forbidden, because energy and momentum con- 
servation cannot simultaneously be fulfilled in this two-body 
process [Fig. 4(c), left panel: the momentum transferred by 
the photon is negligible, but its energy is large]. (This is in 
contrast to high-density electron systems like Gd [8], where 
the momentum necessary for a photon absorption can be pro- 
vided by the multitude of electron-electron scattering events 
occurring during the absorption time [Fig. 4(c), right panel].) 



On the other hand, resonant pumping of the 4/ 7 — > 4/ 6 5c/f2? 
transition is possible. 

Hence, for the first few ps after the photoexcitation, the car- 
rier distribution is comprised of the chemically doped carri- 
ers in the 5d6s conduction band, which remain nearly unaf- 
fected by the pump pulse, and of the fraction of the 4/ elec- 
tron population which is excited into the state with 5d elec- 
tron character. All the mobile carriers contribute to a stronger 
RKKY-like magnetic coupling between the Eu 4/ spins. On 
the other hand, in addition to increasing n c , the Gd doping in- 
duces many-body effects. In particular, low-energy spin fluc- 
tuations in the Gd 5d orbitals enhance the density of states 
at the Fermi level and lead to a downward shift of the con- 
duction band [23], i.e., to a reduction of the 4-f-5d energy gap 
Egap. We have calculated the RKKY-like coupling for the non- 
equilibrium electron distribution described above [24], com- 
puting the doping-dependent many-body renormalization of 
the electronic spectral density as in Ref. [23]. The result is 
shown in Fig. 4(b). Considering the crudeness of the theoreti- 
cal model, the agreement with experimental result in Fig. 4(a) 
is compelling. 

Apart from the semi-quantitative, numerical agreement, the 
non-monotonic behaviors in Figs. 4(a) and 4(b) can be un- 
derstood as follows. The fraction of photoexcited carriers 
is proportional to their spectral density in the final state, 
N(Epf l — Egap), Hence, as E gap is reduced with increasing 
n c (see above), N(E„f, — E gap ) grows. This explains the in- 



10 



°0 -3 

J 'cm J . 



crease of AJ ex {n c )/J ex (n c ) with n c for n c < 10 19 
For n c > 10 20 cm" 3 , AJ ex (n c ) / 4 X ] {n c ) and A/ MSHG // MSHG 
get quenched to negative values for two reasons: (i) With fur- 
ther downward shift of the conduction band the average wave 
number of the carriers increases for fixed excitation energy 
(E p i, = 1.55 eV), so that J ex acquires stronger and stronger 
antiferromagnetic RKKY-like contributions, (ii) With the in- 
creasing metallic nature the electron-electron scattering time 
is reduced, leading to ultrafast thermalization and concomi- 
tant magnetization quenching — conventional demagnetiza- 
tion according to the 3TM begins to dominate. This is nicely 
reflected by the change of the buildup time from 3 ps (man- 
ifestation of RKKY interaction) to ~100 fs (thermalization 
according to the 3TM) in the inset of Fig. 3. 

In summary, we have demonstrated that the ultrafast mag- 
netic coupling dynamics can be tuned from photoinduced en- 
hancement to photoinduced quenching of ferromagnetic order 
in Eui_ r Gd v O by controlling the carrier density. The largest 
enhancement of the ferromagnetic order and the crossover to 
quenching were observed around n c ~ 10 19 and ~ 10 20 cm~ 3 , 
respectively. This behavior is explained by a nonequilibrium 
theory going beyond the established 3TM. Our experimental 
results and their theoretical modelling not only demonstrate 
how to control the stability of a ferromagnetic state at the ul- 
trafast time scale but also give a guide to the material selection 
and design for the dynamical optical control of magnetic or- 
der: Systems with a low density of conduction band carriers 
(which is tuneable) but with a high density of magnetic mo- 



merits (in contrast to diluted ferromagnetic semiconductors) 257401 (2012). 

are in general favorable candidates. 
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